Wave Function Shredding by Sparse Quantum Barriers 
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We discuss a model in which a quantum particle passes through 8 potentials arranged in an in- 
creasingly sparse way. For infinitely many barriers we derive conditions, expressed in terms ergodic 
properties of wave function phases, which ensure that the point and absolutely continuous parts are 
absent leaving a purely singularly continuous spectrum. For a finite number of barriers, the trans- 
mission coefficient shows extreme sensitivity to the particle momentum with fluctuation in many 
different scales. We discuss a potential application of this behavior for erasing the information 
carried by the wave function. 



PACS Nos: 3.65.-W, 2.30.Sa, 5.45.-a, 5.60.Gg 

The prospective use of quantum states as information 
carriers put new emphasis on tailoring systems which 
could process the wave function in a desired way. One of 
such devices, potentially of crucial importance, is an in- 
formation eraser, or briefly shredder. If we look for ways 
in which a quantum system could achieve this goal, a 
natural idea is to seek situations where a slight change of 
the input is magnified significantly in the output. 

In classical mechanics, the sensitive dependence on the 
initial condition of that kind is the essence of the deter- 
ministic randomness exhibited by non-integrable dynam- 
ical systems 0. In quantum mechanics, on the other 
hand, such sensitive dependence has been an elusive ob- 
ject despite many efforts, and a widely accepted coun- 
terpart to the mentioned classical behavior is missing. 
Instead, "quantum signatures" of chaos are often studied 
which typically manifest themselves through statistical 
properties of spectra of various observables There 
has been several attempts to find sensitive dependence 
in quantum mechanics using time periodic models 
However, they are not easy to solve and the numerical 
results often contain ambiguities. In this letter, we aban- 
don such traditional schemes, and instead propose a sim- 
ple one-dimensional model with a potential which has 
many different length scales. 

We start with an infinite number of barriers whose dis- 
tances grow, and ask for conditions under which its spec- 
trum exhibits irregularity, defined in mathematical terms 
as being singularly continuous. We will have unexpected 
encounter with an associated classical system whose er- 
godicity brings about the appearence of purely singularly 
continuous spectrum. We then consider a scattering sys- 
tem in which the number of barriers is finite. Through 
numerical analysis, it is shown that the transmission co- 
efficient exhibits a highly irregular behavior. We also 
show that a wave packet passing through the system un- 
dergoes randomizing during which the shape information 
contained in the input wave function gets lost. 

We choose an explicitly solvable model of a particle 



on the line interacting with (2N + 1) 6 potentials having 
the same coupling constant, which is described by the 
stationary Schrodinger equation 



dx 2 
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ip(x) = k 2 ip(x), 



(1) 



where the positions indexed in ascending order. 

Properties of such systems depend crucially on the po- 
sitioning of the S interactions. In our case they are ar- 
ranged symmetrically, X- n — ~x n with xq = 0, and in 
such a way that their positions are increasingly sparse; 



X\n\+2 



> X 



\n\ + l ~ %\r 



(2) 



We are particularly interested in the situations where the 
distances grow polynomially, exponentially, or faster, 



,1 = b\nf + ce a|n|T 



(3) 



for n 7^ 0, with suitable positive a, b, c, and /3 > 1, 7 > 1. 

We employ transfer matrices which relate solutions of 
Eq. (0) at different points. We define the two component 
vector from the wave function and its derivative as 



ip(x) 



then the transfer matrix M(x,y) is given by 
$(x) = M(x,y)$(y). 



(4) 



(5) 



Because of the zero-range nature of the interaction ||, 
the transfer matrix A4(x,y) is expressed as interlaced 
products of ones describing the free motion with momen- 
tum k, 



M t r(x,y) = 



cos k(x — y) Tsmk(x — y) 



-k sin k(x — y) cos k(x — y) 
and the point interaction transfer matrices 



(6) 
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M(x n + 0,x n - 0) ee M v 



We do not index the latter since the 8 potentials are sup- 
posed to be identical. If x > y the component matrices 
of A4(x, y) are multiplied in the order reversed to that of 
the involved intervals. 

To express the amplitudes for the left-to-right scatter- 
ing over a finite number of barriers, we suppose 



ip(x) = c lkx + r(fc)e 



-ikx 



X < X~N 



t(k)e ikx x > x N 



(8) 



and match the boundary values at x±n using the matrix 
A4 = A4(xn + 0,x_jv — 0); then an easy computation 
using det M. = 1 gives, in particular, 



t(k) 



2ik 



M 2 i - ik(M 2 2 + Mu) - k 2 Mi2 



(9) 



Before making use of the last formula, let us look what 
happens in the limit N — > 00 when the scattering loses 
meaning; we focus on the change of the spectral char- 
acteristics of the system. Traditionally, most attention 
has been paid in physics literature to two types of spec- 
tra the point and the absolutely continuous. Recently, 
there has been a growing recognition of the importance 
of the third type, which is mathematically well studied, 
the singularly continuous spectrum. Here, we ask our- 
selves under which condition the spectrum of our system 
is purely singularly continuous. In view of the symmetry 
it is sufficient to consider the operator 



d 



71=1 



on .£/ 2 (0, 00) with the boundary condition 
ip(0) costf + yj'(0)sintf = 0, 



(10) 



(11) 



in particular, our case with the 6 potential of strength v 
at the origin corresponds to d — — arcsin(l + ^-)~ 1 ^ 2 - 

Several recent papers investigated the relations be- 
tween the spectrum of the operator ( jic| ) and growth prop- 
erties of the corresponding transfer matrices. A sufficient 
condition for the absence of the point spectrum is Q 



dx 



\\M{x,Q)\ 



(12) 



The norm of the transfer matrix can be estimated by the 
product of norms of its component matrices, 



\\M{x ni Q)\\ < v 2 (max(fc, AT 1 ))™ = M[v,k] 



with 
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T + 1 



(13) 



(14) 



Denoting A r 



dx 



x n -\, we have 

N 
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\\M(x,0)\\ 2 ^II^Mlll 2 " 



E 



(15) 



so the integral diverges, for example, if one has A„ > 
cin -1 / 2 e an for some c\ > and 



> 2 {in max(fc, k 1 )+ln^} 



(16) 



If we want a condition independent of the parameters we 
have to require a slightly faster growth, for instance, that 
of (|^) with any 7 > 1. 

Sparse potentials with growing A„ are known to be 
good candidates to produce singularly continuous spec- 
tra. Already the first example by Pearson M is of this 
kind, many others can be found in J^,p. In distinction 
to a typical Pearson potential, however, the height of the 
barriers does not decrease here. This makes it more dif- 
ficult to exclude the absolutely continuous spectrum. To 
show where is the core of the problem we employ the 
modified Prufer technique ||] . It is based on the Ansatz 

tp(x)=R(x)sm9(x), <p'(x) = kR(x) cos 0(x) . (17) 

The amplitude R(x) is constant in the absence of a po- 
tential, so the solution of the Schrodinger equation in the 
interval (x n -i,x n ) between the barriers is 



(p(x) = Rn sin [6 n + k(x - x„-i) 



(18) 



The 6 coupling at the point x n means the wave function 
continuity together with 



(p'(x n + 0) - (p'(x n - 0) = vip(x n ) . 



(19) 



This yields a system of equations for R n ,R n+ \ which is 
solvable under the Eggarter condition PJl0[| 



cot n +i = cot [0 n + kA n ] 



k ' 



(20) 



This determines the phase of the wave function in an 
iterative way. For the ratio of the amplitudes r n := 
R n +i/Rn, we get 



2 _ sin 2 [0 n + kA n ] 

r n ~ -2/1 

sin n+ i 
with the definitions 



= (i + \J [i 2 - 1 sin /3„ 



and 
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(21) 



(22) 



(23) 



The new phase variable (3 n is determined by the evolution 



2 



cot I fcA„+i = cot 



fe' 



(24) 



By Theorem 1.1 of g), the absolutely continuous spec- 
trum will be absent if ||.M(x, 0)|| tends to infinity at least 
for some sequence of points and a.e. k. Since 



||$(x)|| 2 > min(l,fc)i?(Y) 2 , 



(25) 



it is enough to demonstrate the growth for a subsequence 
of {R n }, or equivalently, that the quantity 



N 



(26) 



is exploding for a subsequence of the indices N. Now, for 
a continuous function F S C([— tt, it]), one has 



1 f" 



dft 



(27) 



if the sequence {P n }^Li is uniformly distributed in [—tt, tt] 
jnj. In our case, we have F(f3 n ) := lnr 2 . Using the fact 
that the sine is an odd function we find 



F((i) d/3 



In (fi 2 cos 2 j3 + sin 2 0) d/3 > (28) 



for any \i > 1. It follows that the sequence {i?Ar} is expo- 
nentially growing if {(3 n } is uniformly distributed. Thus 
the problem of singularly continuous quantum spectrum 
is turned into that of dynamical properties of the classical 




FIG. 1. (a) The histogram shows the distribution of the 
phase j3 n in the interval [—tt, tt) for v = k = 1 after N = 20000 
iterations, (b) The distribution of 9 n for the same values. 

We shall not attempt to prove this property and limit 
ourselves to illustrating it for x n — n 2 . The results are 
shown in Fig. 1; similar behavior may be observed for 
other values v, almost all |1^] values of the momentum 
k, as well as for x n — n a with a > 1. Although the 
original phases 9„ may not be uniformly distributed as 
shown in the Fig. 1, it is sufficient that {9 n } is ergodic, 
if the sequence {2fcA„} is uniformly distributed, which 
turns out to be the case when k is an irrational multiple 
of tt, for example, for x n = n 2 [ pd[ . 

Let us sort out the above discussion. If the number 
of barriers is infinite the absolute continuity for Hamilto- 
nians of the type (Q) requires a particular arrangement, 



e.g., a periodic one. On the other hand, randomly dis- 
tributed 6 barriers yield almost surely a pure point spec- 
trum by Anderson localization |l3| . The exclusive char- 
acter of absolutely continuous spectra makes easier to 
get rid of them; the above results suggests that a power- 
like sparseness may be sufficient. At the same time, the 
borderline between dense point and singularly continu- 
ous spectra is extremely unstable and an exponential or 
faster growth may be needed to get a purely singularly 
continuous spectrum. 

Let us return to systems with a finite N. When the 
number of barrriers is finite, the positive energy spectra 
of the system is purely absolutely continuous. However, 
one expects certain trace of singular behavior in observ- 
able quantities even with finite N. One suspects, for 
example, that the poles of the scattering matrix is dis- 
tributed differently for sparse potentials, which should be 
reflected in the increasingly wild behavior of scattering 
amplitudes as one increase the sparseness of the potential 
with fixed finite N. 




FIG. 2. Transmission probability and transmission ampli- 
tude over 2N + 1= 11 barriers for different sparseness. 

That expectation is confirmed in Figure 2, in which we 
plot the transmission coefficient |i(fc) 2 | as a function of 
incident momentum k. The number of barriers is cho- 
sen to be 27V +1 = 11. One has increasing sparseness 
from top to bottom. Irregular fluctuation of the trans- 
mission coefficient which is reminiscent of the classical 
chaotic scattering jy is clearly observed for sparse poten- 
tials. Note also the "ergodic" behavior of the trajectory 
of {Rct(k),lmt(k)} with increasing sparseness. In Figure 
3, we plot the transmission coefficient for a larger num- 



3 



ber of barriers 2iV +1 = 21 with exponential sparseness 
x n = e™. The top graph is magnified by factor of 10 in 
the second graph, part of which is magnified again by fac- 
tor 10 in the third graph. From this, one can discern the 
fact that the irregular fluctuation of quantum scattering 
matrix takes place in many diffrent scale. Notice that 
the despite the irregularity the systems of sparse barriers 
discussed here still exhibit a considerable transmission: 
for the same number of randomly positioned barriers An- 
derson localization takes place and |t(fc)| 2 is practically 
identical to zero. 
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FIG. 3. Transmission amplitude |t(fc)| 2 for exponentially 
sparse 21 barriers in diffrent k scale. 

In actual experimental settings, it is often difficult to 
prepare the particle beam of single wavelength. One 
would rather consider the problem as the transmission 
of wave packet. When one has the incoming wave with 
Fourier component ip(k), the outcoming wave will have 
the spectra t(k)tf(k). The initial wave packet ipi n (x,t) 
will be proccessed after the transmission to give the out- 
comming wave packet 



<fiout(x,t) 



t(x - y)ipin(y,t)dy, 



(29) 



where t{x) is the Fourier transform of t(k). Because 
of the fluctuation in all scale is present in t(k), the 
quantity t(x) has a slow decay at large distances. As 
a consequence, any information contained in the shape 
of the wave packet tpi n (x,t) is effectively washed out in 
fout{x,t) after the convolution (|29|). 

In summary, we have analyzed transport of a quantum 
particle through a sparse family of 6 barriers of identical 
strength. If the barriers are arranged in a non-random 
way there is a non-negligible transmission. When the 



barriers are placed sparsely, the transmission coefficient 
shows an irregular dependence on incoming momentum 
in self-similar manner, and the shape of wave packets 
passing through the system is randomized. Properties of 
the transport can be studied through spectral properties 
of infinite barrier systems. Because of the solvability of 
the problem, a rather detailed analysis can be carried 
out, and a hidden relation between the spectral proper- 
ties of the system and the ergidicity of a certain classical 
dynamical system is revealed. We stress that the simple 
model considered here is not so far removed from ex- 
perimentally realizable system with existing technology. 
Thus we conclude that the realization of initial-condition 
sensitive, information erasing process in purely quantal 
setting may be already within our reach. 
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